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ABSTRACT 

There has been significant recent work by several authors which examines a situation where 
a thin magnetic layer is 'draped' over a core merging into a larger cluster; the same process also 
appears to be at work at a bubble rising from the cluster centre. Such a thin magnetic layer 
could thermally isolate the core from the cluster medium, but only if the same shear process 
which generates the layer does not later disrupt it. On the other hand, if the magnetized layer 
can stabilize against the shear instabilities, then the magnetic layer can have the additional 
dynamical effect of reducing the shear-driven mixing of the core's material during the merger 
process. These arguments could apply equally well to undcrdcnsc cluster bubbles, which would 
be even more prone to disruption. While it is well known that magnetic fields can suppress 
instabilities, it is less clear that a thin layer can suppress instabilities on scales significantly 
larger than its thickness. 

We consider the here stability imparted by a thin magnetized layer. We investigate this 
question in the most favourable case, that of two dimensions, where the magnetic field can most 
strongly affect the stability. We find that in this case such a layer can have a significant stabilizing 
effect even on modes with wavelengths A much larger than the thickness of the layer I - to stabilize 
modes with A « 101 requires only that the Alfven speed in the magnetized layer is comparable to 
or greater than the relevant destabilizing velocity - the shear velocity in the case of pure Kelvin- 
Helmholtz like instability, or a typical buoyancy velocity in the case of pure Rayleigh- Taylor like 
instability. We confirm our calculations with two-dimensional numerical experiments using the 
Athena code. 

Subject headings: hydrodynamics — MHD — instabilities — galaxies: clusters: general — X- 
rays: galaxies: clusters 



INTRODUCTION 



It has been known for some decades in the space science community that an object moving super- 
Alfvenically in a magnetized medi um can very rapidly sweep up a significant magnetic layer which is then 



draped' over the projectile {e.g., iBernikov and Semenovl 119801). After r ecent observations of very sharp 



'cold fronts' in galaxy clusters (see for instance Markevitch and Vikhlinii] ( 2007)), there has been s ignificant 



interest in applying this idea of m agnetic draping in galaxy clusters (e.. 



Vikhlinin et al. 



2001 



2006; 



{e.g.. 



Lvutikov 



Asai et alJl2004ll2005ll2006l) as such a magnetic field could inhibit thermal conduction across the front 



Ettori and Fabian 



2001 



allowing it to remain sharp over dynamically long times. 
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The effect of a strong draped magnetic layer could be even greater for underdense o bjects, such as for 



bubb l es moving through the intercluster medium, as seen in many cool-corc clusters (e.g.. iMcNamara et al 



2005t iBirzan et al.ll2004[). In this cas e, the bubble would be quickly disrupted on rising absent some sort 



of support (e.g.. I Robinson et al.ll2004l ). However, t he draping of a pre-existing magnetic field may strongly 



alter the dynamics, as seen recently in simulations ( Ruszkowski et al. 2007bl la). 



However, the same shear motion which gives rise to the magnetic draping can also drive instabilities 
which could then disr upt the layer. It is clear that magnetic fiel ds can stabilize agains t shear instabilities 
(jChandrasekhaii 1 1 98 11) . and this was considered in this context in lVikhlinin et al.l l|200ir ). However, in that 
work, the instability was considered between two semi-infinite slabs with differing magnetic fields ~ that is, 
the geometric thinness of the magnetic field was not taken into account. This is an appropriate regime for 
considering perturbations much smaller than the thickness of the layer, but less so for modes which could 
disrupt the layer and contribute to the stripping of the core. For these modes, clearly that the layer is in 
fact thin must have some effect, and naively one might expect that the thin layer would only be effective in 
stabilizing modes with size comparable to the breadth of the layer. 

Here we consider the linear growth of Kelvin-Helmholtz and Rayleigh- Taylor instabilities in the presence 
of a thin magnetized layer. We consider the instability in two dimensions, with field lines lying in the plane 
and parallel to the interface, and with perturbations along the direction of the magnetic field. This is the 
case in which the layer could most strongly influence the stability; if such a thin layer could not stabilize 
the flow in such a restricted geometry, it would surely be torn up by instabilities in the more realistic three- 
dimensional case. In three dimensions, interchange modes with wavenumbers perpendicular to the field lines 
are essentially unaffected by the presence of the magnetic field, and thus the field cannot stabilize these 
modes. However, even in this case, overall mixing can be reduced by the presence of a quite weak field {e.g., 
Gardiner and Stond 120071 ) . much weaker than expected in this case, and the introduction of such a stark 
asymmetry (modes in one plane attenuated but in a perpendicular plane unaffected) could have interesting 
observable consequences. 

In [2] we pose the problem, in and S] we derive the growth rates and stability boundaries for the 
Rayleigh- Taylor and Kelvin-Helmholtz instabilities, respectively, a nd in j}5lwe confirm our a nalytic results 
with two-dimensional numerical experiments using the Athena code h.r^^er and Stou . We conclude 



LINEAR THEORY 



We follow the approach and notation of IChandrasekhaij (jl98l[ ). particularly § 105. We begin with the 
equations of two-dimensional, incompressible, inviscid, magnetohydrodynamics: 
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dt 
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where Xi is the i-coordinatc and can take the values (x, z), Ui is the velocity in the i-coordinate direction, p 
is density, 6p is any fluctuation in the density, p is pressure, B is the magnetic field, g is the gravitational 
acceleration with gravity pointing 'down' (e.g., in the direction of — z). and summation over repeated indicies 
is implied. We consider velocity, magnetic fields, and pressure of the form 

U = {u + U,w) u,w <^U, (5) 
B = + b,,h<^B (6) 

where B, W, and p are constant within each region. We will assume that all velocities are highly non- 
relativistic, so that displacement currents and relativistic effects may be neglected. We also assume that the 
flow velocities are much less than the sound speed, and thus may consider incompressible flow. A sketch of 
the situation under consideration is shown in Fig. [TJ 



Following IChandrasekhan (1981) §105, we consider plane wave perturbations of the form e*(*''^+"*) for 
bx, hz, u, w, dp, Sp. Linearizing the equations above results in 

ip [n + kU) u + pU'w - ^B' = -ikdp, (7) 

ip{n + kU)w~^iikbz~b') + 4^pB' = dp' - gSp, (8) 
47r Anp 

ibx{n + kU) ^ ikBu - wB' + bzU' , (9) 

ibz{n + kU) = ikBw, (10) 

iSp{n + kU) ~ —wp , (11) 

iku^w' = 0, (12) 

ikbx + b'^ = 0. (13) 

Because we allow the magnetic field strength to vary between layers, there are extra terms in Eqns. I7I8I9I 
proportional to B' which do not vanish. Within the layers themselves, these terms arc identically zero; but 
their existence will lead to more complex boundary conditions at the interfaces. 



2.1. Solution Within Uniform Layers 

Eauations l9ll0llll andfT^can be solved to express b^jb^^b'^^b'^, Sp, and u in terms of w and its derivatives, 
the system parameters B, B' , p, p', U, the growth rate n, and the wavelength of the disturbance k. This means 
that we can begin to write the dispersion relation n{k) in terms of only the perturbed z-velocity w and the 
known parameters. We can do this by using Eqns. [7] and [H] to eliminate dp, and then eliminate all other 
perturbed variables in favour of w, leaving 

B B' 

2k- (~kwU' + (n + kU) w') + 

47r/9 

in + kUf (—( , + V'kw - (n + kV) w' \ + ((n + kU) (k^w - w") + kwU")^ + 

\ p \(n + kU) ) ^ ' ) 

1.2 o2 

{kw {{n + kU) {nk + U") + /c^C/^ - ikU'"^) - (n + kU) {-2kU'w' + {n + kU) w")) = 0. (14) 



(n + kU) Airp 

Within the three layers, where B' ~ 0, U' = 0, and p' = 0, we have 

,B^ 
Airp 



k^-^-in + kUy ] (k'w~w") =0 (15) 
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Z ->+oo 




Z ^ -00 

Fig. 1. — A sketch of the problem under consideration. We consider a three-layer problem, with two 
semi-infinite, uniform, magnctic-field-free regions separated by a layer with a horizontal magnetic field. We 
consider only field (and disturbances) in the direction of the shear. Wc will generally consider the magnetized 
layer to have the same density and velocity as the top fluid. Disturbances are modeled as periodic in the 
horizontal direction, and vanishing at ±oo in the vertical direction. 
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which, since we are only interested in negative n^, leaves us with only the solutions 

w{z) - e+'=^ w{z) - 6"*=^ 
as in the case of uniform field, and as anticipated in Fig. [1] 



(16) 



2.2. Boundary Conditions Between Layers 



The normal displacement of the interface is -~iw/{n + kU); this can be seen by expressing the perturbed 
interface position Sz in the same form of the other perturbed variables, and the linear-order evolution 
equation for dz becomes, i{n + kU)5z = w. 

Because the displacement of the interface must be unique, w/{n + kU) must be continuous across the 
interfaces. This and a related quantity, the shifted time derivative (n + kU) occurs frequently enough that 
it is useful to express Eq. [T3]in terms of these quantities. Doing so results in 

/ BB' 



47r V " + kU 



+ k^ (n + kU) p 



+ kU 
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n 



kU 



An 

p' = 



(17) 



Integrating this equation over an infinitesimal region across cither of the interfaces gives us the boundary 
conditions across this interface. In doing so, terms that contain no derivatives vanish in the limit, and we 
are left with 



gk' [p] 



n + kU 



p{n + kUy 



n + kU 



B' 



Att \n + kU 







(18) 



where [/] indicates the jump in a quantity / across the interface, and a subscript i refers to a value at 
the interface. Note th at when B is constant, this reduces exactly to the homogeneous field case found in 
Chandrasekiml |l98l[) §106. 



2.3. Matching the Solutions 

In the top and bottom layer, one of the two solution branches (w ^ exp (+fcz) and w ^ exp(— fcz), 
respectively) are clearly unphysical. In the middle layer, however, both can coexist, and so we have as forms 
for the solutions 

w^e~^^ top layer 

W2+e^^^ + 'W2-e~^^ middlelayer (19) 
1/716+'^^ bottom layer 



Because w/{n + kU) must be continuous across the interface, we have 

Wi _ W2+ + W2- 

{n + kU) ~ {n-kU) ' 



(20) 



(n-kU) ^ {n-kU) ■ ^ ' 
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These can be solved for the components of the intermediate velocity in terms of the outer layer velocities 
giving 

wi {n — kU) — W3 (n + kU) 

wie^ ^\n - kU) - W3.{n + kU) 
(g2fci _ i)[n + kU) 



^2+ = 7:^1 TTTZ-TTTv. ' (22) 

^2- = -r^i TTTTTTTT^ • (23) 



n - 



fcC//,,^, ^ "^fe e2.._l „_fct/ e^^'-l 

6^'=' + 1 -71 + fciJ ^ , 1 

We now have two boundary conditions to satisfy — Eq. [TS]at the two interfaces between the layers. 
The top interface condition gives us 

and the bottom interface gives us 



(Ptop - Pbot) (^^^ - (^Ptop (n - fcC/) (;^) ^^^^^ - Pbot^^ifc j (^^^ j ^^^^ = (27) 

This gives us two equations in terms of w\ and 1^3. Using our expressions for wi± and then solving the 
first equation for W3 in terms of and substituting the result into the second gives us our final dispersion 
relation of 

({n - kUf - i/fc2„2 A ^ ^2j^2 ^ Atfc {g - 2nU) - ^Pvl) 
(e^''^ - 1) ^ , = (28) 

where the Atwood number, At, is defined to be the non-dimensional density difference 

At = ^^°P ~ ^'^"^ . (29) 

Ptop + Pbot 

It is worth noting that in the no magnetization limit va 0, the only term containing I can be divided 
out, so that the solution does not depend on I; this is as it must be, as without a magnetic field nothing 
distinguishes the ^-thick middle layer from the top layer. Also, in the infinitely thin layer limit exp 2kl —>■ 1 
and the solution reduces to the non-magnetized case. 



3. EFFECT ON THE RAYLEIGH-TAYLOR INSTABILITY 

If we consider the 'pure' Rayleigh- Taylor instability, with no horizontal shear, then U ^ and we are 
left with 



„2 + Atg/fc + (e^'^' - 1) V" -^fcM)(n^ + Atgfc~lj^fcM; ^ ^ ^g^^ 

7^2 — fc^V^ 
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This equation relates the growth rate n to two other inverse timescales on the scale of - an Alfven 
frequency = k'^v\, and a gravitational timescale t^'^ — Kigk. Expressing the growth rate and Alfven 
frequency in units of the inverse gravitational timescale, and ignoring the trivial Alfven wave solution n = u>a, 
the result is a quadratic in n^: 

g2fei^4 I I l-i , ,2 \ I /'„2A:; ^\ ( ^ ^^"^^,.2\\~2 , ,2 ( ^ i ^ ( ^2kl --^l. 1 + At 



+ \{X + (e^^-' ^-ij j n'-^l [I + ^ {e'"' ^ l) (l - ^-1) ) = 0. (31) 

To consider the degree to which the magnetized layer stabilizes against the Rayleigh- Taylor instability, we 
consider (for simplicity) the maximally unstable case where At ^ 1. For stability, it is necessary and 
sufhcient that the roots of the quadratic in be positive and real; a quadratic of the form ax^ + bx + c ^ 
has positive and real roots for c/a > 0, b/a < 0, and — Aac > 0. The first two of these conditions reduce 
to 

2kl I -1 

> ^ (32) 

2 2e2fc' , ^ 

'^A > 3g2M „ 1 ' (33) 

and the third is satisfied for all real u!\. Of these conditions, the first controls, as the second simply posts 
a lower limit for uj\ of between 2/3 and 1, while the lower limit for the first is always greater than 1. Thus 
the condition for stability in the case of a Rayleigh- Taylor instability with At = 1, which would otherwise 
always be unstable, is 

-A > (5^) 9k-' (34) 

with no strength of magnetic field able to stabilize in the case of an infinitely thin layer. This stability 
criterion is shown in Fig [21 



4. EFFECT ON THE KELVIN-HELMHOLTZ INSTABILITY 

The shear terms in the case of the Kelvin- Helmholtz instability make the dispersion relation significantly 
more complicated, 

n"* - 2 (1 + At) kUn^ + ie-^'^' (Atk^v\ + e^^^ (4 (1 + 2At) k^U^ - (2 + ki)k'^v\)) r?-V 

-e-'^^^kU {k^v\ + e^^^ (2(1 + ki)k^U'^ - (1 2k\)k^Vj)) ri+ 

ie-2*^-'(2Atfc2[/2fc2^2 _(^^A^)^4^4 ^g2fei(2fc2c/2_fc2^2)(2fc2c/2_(i + At)^^^ - 0. (35) 

For considering the stability boundary we will again consider the most unstable case, in this case At = 0. 
As in the previous case, the growth rate and the Alfven frequency can be expressed in terms of the other 
timescale of the problem, the advection time across the wavelength of the perturbation, r^/ = 1/ (kU), leaving 
us with 

n"* - 2h^ - {uj\ -2)n^ + {u;\ (l - e^^fci) _ 2) n 

+ \[{u.\-2)'~c.\e'''^) = 0. (36) 
While this is a fairly unpleasant quartic in n, it is a fairly approachable quadratic in w^: 

i (e^^' - 1) + (1 + n (e-2'^-' - l) + n^) uj\-{fi- lf{f^ + 1) = 0. (37) 
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Fig. 2. — Plotted is the Rayleigh- Taylor stability boundary for the most-unstable (At = 1) case, as given 
in Eq. 1341 The magnetic field strength in the magnetized layer (here expressed as the Alfven speed squared 
in units of g/k) necessary to stabilize against a mode is plotted as a function of the thickness of the layer 
in units of the wavelength of the mode. Thus, to stabilize modes an order of magnitude greater longer than 
the layer is thick (//A ~ 0.1), a magnetic field such that v\ > 2g/k is necessary. 
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The two roots are 



- 2 g2H „ 1 ■ (38) 

For stability, wc consider the neighborhood around Im{fi) = 0. In this case, both these branches have 
minima for purely oscillatory modes around n = 1; but in this neighborhood the second, negative, branch, 
has no real solutions for uja with Im(n) ^ 0, so cannot be relevant to the question of stability. The positive 
branch has a minimum of 

2kl I 1 

-A = 4^5m^ (39) 
so that stability is ensured when this condition is met, or 

2 . ,'e2^' + l 



< > [-^^^ ) i2Ur, (40) 

which is the same condition for stability of the Rayleigh- Taylor instability, but with 2U replacing g/k; recall, 
however, that the two conditions arc for two different values of the Atwood number, such that the instability 
is maximized in each case. This stability criterion is plotted in Fig. [31 



NUMERICAL RESULTS 



To confirm the results of the pre vious sections, numerical e xperiments were performed in two dimensions 
using version 3.0 of the Athena code ( Gardiner and Stonel2005 ). a dimcnsionally imsplit, highly configurable 



MHD code. For the results in this work we used the ideal gas MHD solver with an adiabatic equation of 
state (7 = 1.4), and the 3rd-order accurate solver using a Roe- type flux function. We considered a domain of 
size, in code units, of [—1/6, 1/6] x [—1/6, 1/6] with resolution 400 x 400 for the Kelvin-Hclmholtz instability 
simulations, and with a slightly vertically extended domain ([—1/6,1/6] x [—1/4,1/4], 400 x 600) for the 
Rayleigh- Taylor instability experiments. 

To ensure that the resolution used was adequate, a resolution study was performed on a fiducial run (a 
Rayleigh- Taylor simulation with B^fi = 0.07, va = 0.0495; the magnetic field code units are such that the 
Alfven speed, v\ = i?^/2) with resolution varying between a factor of two less than this resolution and a 
factor of two more; measured growth rates varied by only approximately ±3 percent. 

The analytic results presented in previous sections were in the incompressible limit; in our numerical 
experiments here the fiducial density was 1 in code units and the pressure was set so that the sound speed Cs 
would be 1 in code units, which is an order of magnitude larger than the velocities achieved in either set of 
simulations; thus the Mach number M < 0.1, and incompressibility remains a reasonable approximation. In 
both set of simulations, a magnetized layer of thickness 1/60 was initialized starting at y — 0, with strength 
of magnetic field varied from run to run. To keep the initial conditions in pressure equilibrium, the thermal 
pressure was reduced in this layer, but because of the large sound speed (and consequently large plasma (3) 
this was a small reduction (never more than a few percent). 

In both sets of simulations, the interface was given a sinusoidal velocity with a wavelength equal to 
the size of the box, and an amplitude of Vpert ~ 0.0025. In the Kelvin-Helmholtz case, a background shear 
velocity of —0.1 in the x direction was applied to the top layer and the magnetized layer, and of -1-0.1 in 
the bottom layer. In the Rayleigh- Taylor case, a gravitational acceleration of g = 0.1 was applied in the 
negative y direction. Both set of simulations used an Atwood number At = 1/3, so that the top layer had 
density pt = 2. Snapshots of the evolution of the simulations are shown in Figs|3]and[ni 
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Fig. 3. — Plotted is the Kelvin-Helmholtz stability boundary for the most-unstable (At = 1) case, as given 
in Eq. 1401 The magnetic field strength in the magnetized layer (here expressed as the Alfvcn speed squared 
in units of the half-shear velocity squared, U'^) necessary to stabilize against a mode is plotted as a function 
of the thickness of the layer in units of the wavelength of the mode. Thus, to stabilize modes an order of 
magnitude greater longer than the layer is thick (//A ~ 0.1), a magnetic field such that v'j^ > 2(2[/)^ is 
necessary. 
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Fig. 4. — Snapshots from representative runs of the Raylciglr- Taylor instability with a magnetized layer, 
performed with the Athena code. Shading represents density, and contours indicate magnetic field strength. 
In code units, the density of the top fluid is pt = 2, and that of the bottom is = 1 (At = 1/3), and pressure 
is set such that in the bottom layer the adiabatic sound speed is 10, and the domain is in total pressure 
equilibrium. The wavelength of the initial perturbation is 1/3, and the thickness of the magnetized layer is 
1/60; the acceleration due to gravity g is 1/10 and in the negative y direction. Shown are runs with varying 
initial horizontal magnetic fields, with, from top to bottom, B^^q ~ 1/200, 1/20, 1/10. Snapshots are shown 
at times in code units of, left to right, approximately (1.75,3.5,5.25). The contours are for magnetic field 
strength of i?x,o/4 and 3Bx.o/4:. 
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Fig. 5. — Snapshots from representative runs of the Kelvin-Helmholtz instability with a magnetized layer, 
performed with the Athena code. Shading represents density, and contours indicate magnetic field strength. 
In code units, the density of the top fluid is pt = 2, and that of the bottom is pf = 1 (At = 1/3), and 
pressure is set such that in the bottom layer the adiabatic sound speed Cs is 1, and the domain is in total 
pressure equilibrium. The wavelength of the initial perturbation is 1/3, and the thickness of the magnetized 
layer is 1/60; the acceleration due to gravity g is 1/10 and in the negative y direction. Shown are runs 
with varying initial horizontal magnetic fields, with, from top to bottom, Bx.o = 0.1,0.2,0.3. Snapshots are 
shown at times in code units of, left to right, approximately (0.75, 1.5,2.25). The contours are for magnetic 
field strength of i?x,o/4 and iB^.o/^- 
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From the outputs of the simulations, growth rates for the instabihties were calculated by considering 
the growth of the amplitude (measured by finding the mean vertical position of the magnetized layer and 
fitting to a sinusoid of the wavelength of the perturbed mode) as a function of time. Exponentials were fit to 
this scries of amplitudes, for those times where the amplitude was resolved by at least three zones and where 
the amplitude was less than 1/15 of the wavelength {e.g., before nonlinear evolution begins to matter). For 
the simulations reported here, this means the fit was performed with amplitudes in the range [1/400, 1/45], 
covering approximately a decade in amplitude. Plotted in Fig. [6] are four examples of this procedure for the 
Rayleigh- Taylor simulations. For comparison, the results of a resolution study of a fiducial Rayleigh- Taylor 
case is shown in Fig. \7\ 

Because of the low speeds of these flows, many timesteps (typically on order 20000) must be take to 
evolve these instabilities. Over that period of time modest amounts of numerical diffusion slightly modify 
the structure of the magnetic layer, as does the still slightly compressible flow itself; the change is shown in 
FiglH The growth rate of the instabilities is very sensitive to the thickness of the layer, and this modification 
of the magnetic layers thickness must be taken into account when comparing with analytic results. 

The comparison of measured growth rates to analytic results is shown in Figs. [91 and [TUl with theoretical 
results taken from the numerical solution of Eqs. l30l and l35l Where the numerical experiments can give clean 
growth rate measurements, the analytic and numerical results agree to within a few percent, and there is 
also agreement at the few percent level on the boundary of stability. 



6. DISCUSSION 



We have shown, through derivation of dispersion relations and two-dimensional numerical experiments, 
that it is possible for even thin magnetized layers to suppress instability growth on scales much larger than 
their own thickness as long as the magnetic field strength is high enough that the Alfven velocity in the 
layer in the direction of the perturbation is of order the relevant destabilizing velocity scales. In the case 
of the Kelvin-Helmholtz instability, the most most relevant case for core mergers, a magnetized layer can 
stabilize modes an order of magnitude larger than the thickness of the layer if the Alfven speed of the same 
magni tude as the full shear velocity va ~ '2U. But this is almost automatically true; as shown in iLvutikov 
(|2006[ ). near the stagnation line the magnetic pressure reaches equipartition with the ram pressure, meaning 
that this condition on the velocities is met. Thus we would expect, certainly near the stagnation line, that the 
magnetic draping 'protects' a merging core f rom instabi lities, as would appear to be the case in simulations 
published in the literature |Asai et al.ll2004l . boOSl |2006| . Dursi and Pfrommer, in preparation). 



Clearly, this stabilization only applies to perturbations along the field; but in this plane, then, the 
draped field can then have the twin effects of protecting the merging core against thermal disruption, and 
reducing the shear effects which would tend to mix in the core material earlier. In the plane perpendicular 
to the magnetic field, shear-driven instability will occur unimpeded, leading to a distinct asymmetry in the 
resulting magnetic layer and, presumably, the moving core or bubble. Such asymmetries have been seen in 
3d simulations in the literature, as cited above, and would potentially be observable. It is interesting to 
note however that even a quit e weak field can effect global mixing in the presence of similar instabilities 
(e.g.. iGardiner and Stondl2007^ ■ and so even in the plane perpendicular to the field it is possible that a thin 
magnetized layer could keep an interface sharper than could exist absent the magnetic field. 

The full three dimensional stability problem remains to be tackled, and less clear still is the effect of a 
more realistic magnetic field, not expected to be planar or uniform, and the effects of fully three-dimensional 
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Fig. 6. — Evolution of amplitude of the Rayleigli- Taylor perturbation with time, for several runs (plotted 
are, top to bottom, Bq = (0.005,0.05,0.1,0.2)), shown as symbols; also plotted here are best-fit exponentials 
to measure the growth rate of the instability. Exponentials were fit to data from those times where the 
amplitude was resolved by at least three zones and where the amplitude was less than 1/15 of the wavelength 
{e.g., before nonlinear evolution begins to matter); for the simulations reported here, this means the fit was 
performed with amplitudes in the range [1/400, 1/45], covering approximately a decade in amplitude. The 
range considered is indicated by thick shaded lines; below the bottom line the 'jumpiness' due to poor 
resolution on the grid is visible, and above deviations from exponential behaviour due to nonlinear effects 
becomes evident. 
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Fig. 7. — As in Figure [6l but with Bq = 0.07, and varying the resolution. Shown are for several runs with 
varying by a factor of 4, where is the horizontal resolution in the simulation, with vertical resolution 
kept proportionate; thus N is the number of points used to represent one wavelength, with the thickness of 
the layer being represented by 7V/20 points. A^ = 100 is not quite enough resolution to accurately measure 
the growth rate, but the difference in fitted exponential growth rates between 200, 400, 600, and 800 is only 
a few percent. 
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Fig. 8. — Magnetic field profile through a; = for the B^.o — 1/10 Rayleigh- Taylor simulation at time (filled 
squares) and 3.75 (open circles). The later-time profile is shifted leftwards by Sy — 0.00675 for comparison. 
During the course of the simulation, the thickness of the magnetic field layer at its peak is reduced to about 
18 grid points from the original 20. Besides the inevitable effects of numerical diffusion slightly smearing 
out the profile over approximately 18000 timesteps, effects from the slight compressibility of the flow, and 
asymmetry from the direction of gravity (to the left in this plot) is shown by the slight increase in magnetic 
field density and the slope at late times. 











0.02 



0.04 



0.06 

va 



0.08 



0.1 



0.12 



Fig. 9. — Rayleigh- Taylor growth rates as a function of magnetic field strength (given here in terms of Alfvcn 
speed, i'a) in the magnetized layer, measured from simulation (diamonds) and predicted from analytic theory 
(lines). The dashed line gives the prediction using the nominal width of the magnetized field layer, and the 
solid line gives the prediction for the actual width of the magnetized layer in the simulations after diffusion 
eats away at the profile as shown in Fig. [8l A gray point indicates a run done at twice the resolution to 
check to see if the default resolution was adequate. The theoretical curves are taken from numerical solution 
of the dispersion relation, Eq. [30l 
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Fig. 10. — Kelvin-Hclmholtz growth rates as a function of magnetic field strength (given here in terms of 
Alfvcn speed, va) in the magnetized layer, measured from simulation (diamonds) and predicted from analytic 
theory (lines). The dashed line gives the prediction using the nominal width of the magnetized field layer, 
and the solid line gives the prediction for the actual width of the magnetized layer in the simulations after 
diffusion eats away at the profile as shown in Fig. [S] Below va ~ 0.1, smaller- wavelength modes, seeded by 
numerical grid noise, grow very quickly and contaminate any attempt to measure a linear growth rate for 
the mode under consideration. The theoretical curves are taken from numerical solution of the dispersion 
relation, Eq. [35l 



- 19 - 



perturbations on such a layer. Consideration of this more comphcated and reahstic case is left to future 
work. 
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